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We study the phase transition of the three-dimensional complex \ theory by considering the 
geometrically defined vortex-loop network as well as the magnetic properties of the system in 
the vicinity of the critical point. Using high-precision Monte Carlo techniques we examine an 
alternative formulation of the geometrical excitations in relation to the global (9(2)-symmetry 
breaking, and check if both of them exhibit the same critical behavior leading to the same critical 
exponents and therefore to a consistent description of the phase transition. Different percolation 
observables are taken into account and compared with each other. We find that different defini- 
tions of constructing the vortex-loop network lead to different results in the thermodynamic limit, 
and the percolation thresholds do not coincide with the thermodynamic phase transition point. 
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Three-dimensional, globally 0(2) symmetric spin models and field theories exhibit line-like 
topological excitations which form closed networks. An issue of central importance is the question 
whether the distribution of these vortex lines and their percolation properties do indeed encode 
the critical exponents of the thermodynamically defined phase transition. More specifically, the 
question we want to address in this paper is: Is there a similar clue in the case of vortex networks 
as for spin clusters, or do they display different features? Percolational studies of spin clusters 
showed that the geometric approach only works, if one uses a proper stochastic (Fortuin-Kasteleyn) 
definition of clusters |2|, |3l Hj. When connecting the vortex line elements to closed loops, which 
similar to spin clusters are geometrically defined objects, and a branching point with n > 2 junctions 
is encountered, a decision on how to continue has to be made. This step involves a certain ambiguity 
and gives room for a stochastic definition. In particular we want to investigate the influence of the 
probability of treating such a branching point as a knot. 

In this paper we concentrate on the three-dimensional (3D) complex Ginzburg-Landau model, 
the field theoretical representative of the 0(2) universality class, with normalized lattice Hamilto- 
nian 



N -ff„ „ 1 d 



H[Y] = £ [^(|Vn| 2 -l) 2 +2 L \Yn-Yn+fl\ 2 



(1) 



where \\f n = \j/ njc + i\jf ny = \ y n \e"^" is a complex field and a is a temperature independent parameter. 
/I denotes the unit vectors along the d coordinate axes, N = L d is the total number of sites, and an 
unimportant constant term has been dropped. The partition function reads 



jD\i/Dye- H/T , (2) 



where / Di/aDi// = / DRei/fDImi/A stands short for integrating over all possible complex field con- 
figurations. In the limit of a large parameter a, it is easy to read off from Eq. (|l|) that the modulus 
of the field is squeezed onto unity such that the XY model limit is approached with its well-known 
continuous phase transition in 3D at T c = l/j3 c - ~ 2.2 

In order to characterize the transition we performed Monte Carlo simulations and measured 
among other quantities the energy (H), the specific heat c v = ({H 2 } — (H) 2 )/N, the mean-square 
amplitude {\y\ 2 ) = (l/A^)L^=i(IV / »| 2 )' an d the magnetization (M) = Ym=i(Yn)- The main focus 
of this paper is on the properties of the geometrically defined vortex-loop network. The standard 
procedure to calculate the vorticity on each plaquette is by considering the quantity 

m = ~ fahx + [02 - fohn + [03 - 04]2^ + [04 ~ Qlhit) , 0) 

2K 

where 0i , . . . , 04 are the phases at the corners of a plaquette labelled, say, according to the right- 
hand rule, and [a]27r stands for a modulo 2n: [oC\2n = cc + 2%n, with n an integer such that a + 
2nn € (—71, 7t], hence m = nn + «23 +«34 +«4i- If m ^ 0, there exists a topological charge which 
is assigned to the object dual to the given plaquette, i.e., the (oriented) line elements */ jU which 
combine to form closed networks ("vortex loops"). With this definition, the vortex "currents" 
can take three values: 0, ±1 (the values ±2 have a negligible probability and higher values are 
impossible). The quantity v = ^L„,^ |*^u,n| serves as a measure of the vortex-line density. 

In order to study percolation observables we connect the obtained vortex line elements to 
closed loops. Following a single line, there is evidently no difficulty, but when a branching point, 
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where n > 2 junctions are encountered, is reached, a decision on how to continue has to be made. 
If we connect all in- and outgoing line elements, knots will be formed. Another choice is to join 
only one incoming with one outgoing line element, with the outgoing direction chosen randomly. 
We will employ two "connectivity" definitions here: 

• "Maximal" rule: At all branching points, we connect all line elements, such that the maximal 
loop length is achieved. That means each branching point is treated as a knot. 

• "Stochastic" rule: At a branching point where n > 2 junctions are encountered, we draw a 
uniformly distributed random number G (0, 1] and if this number is smaller than the con- 
nectivity parameter c we identify this branching point as a knot of the loop, i.e., only with 
probability < c < 1 a branching point is treated as a knot. In this way we can systematically 
interpolate between the maximal rule for c = 1 and the case c = 0, which corresponds to the 
procedure most commonly followed in the literature [jg]. 

We can thus extract from each lattice configuration a set of vortex loops, which have been 
glued together by one of the connectivity definitions above. For each loop in the network, we 
measure, among others [f7|], the following observables: 

• "Extent" of a vortex loop in 1, 2, or 3 dimensions, ^iu, and ff^: This means simply 
to project the loop onto the three axes and record whether the projection covers the whole 
axis, or to be more concrete, whether one finds a vortex-line element of the loop in all planes 
perpendicular to the eyed axis. This quantity can thus be interpreted as percolation probabil- 
ity [JSj] which is a convenient quantity for locating the percolation threshold j8 p . 

• "Susceptibilities", %i- F° r the vortex-line density v and any of the observables &i defined 
above, one can use its variance to define the associated susceptibility, %i = N({6>f) ~ (^') 2 )» 
which is expected to signal critical fluctuations. 

To update the direction of the field []|], we employed the single-cluster algorithm Jl0| ] similar to 
simulations of the XY spin model P|. The modulus of y was updated with a Metropolis algorithm. 



Here some care was necessary to treat the measure in (Q) properly (see Ref. JTlj]). One sweep 
consisted of N spin flips with the Metropolis algorithm and N sc single-cluster updates. For all 
simulations the number of cluster updates was chosen roughly proportional to the linear lattice 
size, N sc ~ L, a standard choice for 3D systems as suggested by a simple finite-size scaling (FSS) 
argument. We performed simulations for lattices with linear lattice size ranging from L = 6 up 
to L = 40, subject to periodic boundary conditions. After an initial equilibration time of 20000 
sweeps we took about 100000 measurements, with ten sweeps between the measurements. All 
error bars are computed with the Jackknife method. 

In order to be able to compare standard, thermodynamically obtained results (working directly 
with the original field variables) with the percolative treatment of the geometrically defined vortex- 



loop networks considered here, we used the same value for the parameter o = 1.5 as in Ref. [ |12| ] 
for which we determined by means of standard FSS analyses of the magnetic susceptibility and 
various (logarithmic) derivatives of the magnetization a critical coupling of 

& = 0.78008(4) . (4) 
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Table 1: The critical exponents of the 3D XY model universality class as reported in Ref. jl3| ] and the 
correction-to-scaling exponent to of Ref. 



a 


P 


7 


8 


TJ V ft) 


-0.0146(8) 


0.3485(2) 


1.3177(5) 


4.780(2) 


0.0380(4) 0.67155(27) 0.79(2) 



Focussing here on the vortex loops, we performed new simulations at this thermodynamically 
determined critical value, j8 = 0.78008, as well as additional simulations at /3 = 0.79, 0.80, and 
0.81. The latter j8 values were necessary because of the spreading of the pseudo-critical points of 
the vortex loop related quantities. As previously we recorded the time series of H, M, and 
|y/"| 2 , as well as the helicity modulus and the vorticity v. In the present simulations, however, 
we saved in addition also the field configurations in each measurement. This enabled us to perform 
the time-consuming analyses of the vortex-loop networks after finishing the simulations and thus 
to systematically vary the connectivity parameter c of the knots. 

The FSS ansatz for the pseudo-critical inverse temperatures /3,-(L), defined as the points where 
the various %i w& maximal, is taken as usual as 

ft(L) = A-c + dL-Vv + C2L -i/v- ffl + ... f (5) 

where (5j c denotes the infinite-volume limit, and v and ft) are the correlation length and confluent 
correction critical exponents, respectively. Here we have deliberately retained the subscript i on 

A> 

Let us start with the susceptibility % v of the vortex-line density, which plays a special role 
in that it is locally defined, i.e., does not require a decomposition into individual vortex loops. 
Assuming XY model values for v and ft) (cf. Table [T]) and fitting only the coefficients j8,- c and c/, 
we arrive at the estimate j8 vc - = 0.7797(14) with a goodness-of-fit parameter Q = 0.20. This value 
is perfectly consistent with the previously obtained "thermodynamic" result (Q). On the basis of 
this result one would indeed conclude that the phase transition in the 3D complex Ginzburg-Landau 
field theory can be explained in terms of vortex-line proliferation [[Tj]. 

To develop a purely geometric picture of the mechanism governing this transition, however, 
one should be more ambitious and also consider the various quantities €f[ introduced above that 
focus on the percolative properties of the vortex-loop network. As an example for the various 
susceptibilities considered, we show in Fig. |l|(a) the susceptibility X3D of <^ 3D for c = and c = 1. 
The scaling behavior of the maxima locations j83d (L) of the susceptibility X3D of ^?3d for c = and 
c = 1 is depicted in Fig. |](b), where the lines indicate fits according to Eq. ^ with exponents fixed 
again according to Table |. We obtain p 3D;C = 0.7824(1) with £ 2 /dof = 1.14 (Q = 0.32, L > 8) 
for c = and /3 3DjC = 0.8042(4) with £ 2 /dof = 0.75 (Q = 0.58, L > 20) for c = 1. While for the 
"stochastic" rule with c = the infinite- volume limit of ^3d(L) is at least close to j8 c , it is clearly 
significantly larger than /3 C for the fully knotted vortex networks with c = 1. 

By repeating the fits for all vortex-network observables and the parameter c between and 
1 in steps of 0.1, we observed that the location of the infinite- volume limit j3,- c does depend on 
the connectivity parameter c used in constructing the vortex loops in a statistically significant way. 
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Figure 1: (a) Susceptibility of as a function of j3 for the "stochastic" (c = 0) and the "maximal" 
(c = 1) rule, (b) Location of the percolation thresholds determined from the maximum of ^30 as a function 
of L~ l l v . The lines indicate fits using the FSS ansatz (^) with v and CO fixed according to Table |l|. The 
horizontal dashed line shows the thermodynamically determined critical coupling j3 £ = 0.78008(4). 



With decreasing c, the infinite-volume extrapolations come closer towards the thermodynamical 
critical value (Q), but even for c = they clearly do not coincide. 

We nevertheless performed tests whether at least for c = the critical behavior of the vortex- 
loop network may consistently be described by the 3D XY model universality class. As an example 
for a quantity that is a priori expected to behave as a percolation probability we picked the quantity 
^3D- As is demonstrated in Fig. ||(a) for the case c = 0, by plotting the raw data of as a function 
of p for the various lattice sizes, one obtains a clear crossing point so that the interpretation of 0^, 
as percolation probability is nicely confirmed. To test the scaling behavior we rescaled the raw 
data in the FSS master plot shown in Fig. |2](b), where the critical exponent v has the XY model 
value given in Tableland /3 c (^3d) = 0.7842 was independently determined by optimizing the data 
collapse, i.e., virtually this is the location of the crossing point in Fig. ||(a). The collapse turns out 
to be quite sharp. For c > we found also a sharp data collapse, but for a monotonically increasing 
exponent v, which is for large c values compatible with the critical exponent v = 0.8765(16) of 



3D percolation [15]. One should keep in mind, however, that neither p3D lC as extrapolated from 
the susceptibility peaks nor the estimate obtained from the crossing point in Fig. ||(a) is compatible 
with j3 c . 

To summarize, we have found for the 3D complex Ginzburg-Landau field theory that the geo- 
metrically defined percolation transition of the vortex-loop network is close to the thermodynamic 
phase transition, but does not quite coincide with it for any observable we have considered. Our 
results for the connectivity parameter c G [0, 1] extend the claim of Ref. |Q] for the 3D XY spin 
model that neither the "maximal" (c = 1) nor the "stochastic" rule (c = 0) used for constructing 
macroscopic vortex loops does reflect the properties of the true phase transition in a strict sense. A 
more detailed presentation of these and additionally results for several other observables is given 
in Ref. 
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Figure 2: (a) as function of j3 for c = 0. (b) Rescaled data with v fixed at the 3D XY model value 
(cf. Table [l]) and choosing /3 c (^3d) = 0.7842 from the location of the crossing point in (a) for the best data 
collapse. 
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